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ABSTRACT 

The aim of this paper is to investigate spectral and photometric properties of 854 faint (Jab ^ 25 mag) star-forming galaxies (SFGs) 
at 2 < z < 2.5 using the VIMOS Ultra-Deep Survey (VUDS) spectroscopic data and deep multi-wavelength photometric data in three 
extensively studied extragalactic fields (ECDFS, VVDS, COSMOS). These SFGs were targeted for spectroscopy based on their 
photometric redshifts. The VUDS spectra are used to measure the UV spectral slopes (j8) as well as Lya equivalent widths (EW). 
On average, the spectroscopically measured /? (-1.36±0.02), is comparable to the photometrically measured p (-1.32±0.02), and has 
smaller measurement uncertainties. The positive correlation of p with the Spectral Energy Distribution (SED)-based measurement 
of dust extinction Es(B-V) emphasizes the importance of p as an alternative dust indicator at high redshifts. To make a proper 
comparison, we divide these SFGs into three subgroups based on their rest-frame Lya EW: SEGs with no Lya emission (SFGn; 
EW < OA), SFGs with Lya emission (SFGl; EW > OA), and Lya emitters (LAEs; EW > 20A). The fraction of LAEs at these redshifts 
is ~-10%, which is consistent with previous observations. We compared best-fit SED-estimated stellar parameters of the SEGn, SFGl 
and LAE samples. Eor the luminosities probed here (~L*), we find that galaxies with and without Lyo- in emission have small but 
significant dilferences in their SED-based properties. We find that LAEs have less dust, and lower star-formation rates (SFR) compared 
to non-LAEs. We also find that LAEs are less massive compared to non-LAEs, though the difference is smaller and less significant 
compared to the SER and Es(B-V). When we divide the LAEs based on their Spitzer/TRAC 3.6yum fluxes, we find that the fraction of 
IRAC-detected (m 3 5 < 25 mag) LAEs is much higher than the fraction of IRAC-detected narrow band (NB)-selected LAEs at z- 2-3. 
This could imply that UV-selected LAEs host a more evolved stellar population, which represents a later stage of galaxy evolution 
compared to NB-selected LAEs. 

Key words. Galaxies: evolution - Galaxies: formation - Galaxies:high redshift - Cosmology: observations 


* Based on data obtained with the European Southern Observatory 

Very Large Telescope, Paranal, Chile, under Large Program 185.A- 


0791. 


1. Introduction 

In the last decade, large numbers of star-forming galax¬ 
ies (SFGs) — selected using the Lyman break color tech- 
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nique and/or the photometric redshift technique — were 
identified an d studied at z > 3 using deep photometric sur- 
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most of these galaxies are very challenging because of their faint 
magnitudes and lack of high resolution multi-wavelength data. 
Therefore, detailed studies of SFGs at lower redshifts (z < 3) us¬ 
ing ground-based spectroscopy and multi-wavelength photom¬ 
etry are vital for understanding physical processes at the peak 
epoch of star formation as well as to better understand high red- 
shift galaxy formation. 

The study of SFGs at z~2 has a significant impact on 
our understanding of galaxy properties. First, these SFGs are 
at redshifts corresponding to the peak epo ch (z~ 1-3) of the 
global star formation rate (SFR) density (e.g..lMadau etalJ 19981 : 


ICucciati etalJl201^lMadau & Dickinsonll2014 ). where >50% of 

the stars in the present-day Universe formed, and is a perfect 
epoch to study star formation processes. Second, this cosmo- 
logically interesting redshift is also where we have access to a 
wealth of multi-wavelength data, and rest-frame optical spec¬ 
troscopy, to study galaxy properties that we cannot investigate at 
higher redshifts with current technology. The major advantage of 
identifying and studying various physical properties — includ¬ 
ing star formation properties — of these SFGs is that they can 
be investigated in rest-frame UV as well as rest-frame optical fil¬ 
ters. Third, these galaxies are likely lower redshift analogs of the 
high redshift SFGs, and putting in a context such analysis will 
help to she d light on the process of reionization in the early Uni¬ 
verse (e.g.. lLabbe etaklUOlOHStark etal.ll2010h and how the first 
galaxies formed. Finally, the redshift z ~ 2 is the lowest redshift 
where we can identify and study Lya in emission from ground- 
based spectroscopy and these galaxies are bright enough for ex¬ 
tensive spectroscopic studies on 8-lOm class telescopes. 

Spectroscopy is an essential and powerful tool to fully un¬ 
derstand galaxy properties but the underlying photometric se¬ 
lection of galaxies can bias the spectroscopic sample proper¬ 
ties. The primary techniques to select SFG s by color a t z-2 
are: (1) sBzK (using the B, z, K bands, [Daddi et ^ |2004 
2007 ), (2) BX/BM (using t he U, G, R bands, Steidel et^ 


2004 I Adelberger et al.l l2004l) . and (3) LEG (using the bands 


which bracket th e redshifted Lyman limit, iHathi et ^ 120101 : 
lOesch et al.l2010h . The other main approach is based on the mag¬ 
nitude/flux limit and/or photometric redshift selection, as used 
by many spe ctroscopic surveys such as, WPS, GMASS, zCOS- 


MOS fe.g.. iT.e Fevre etal.l 1200.51 iT.illv etlTI 12001 iKiirk etal.1 


l2009t iLe Fevre et al. 20131) . On the other hand, SFGs are also 


selected bas ed on the emissio n -line/Narrow-band (NB) tech¬ 
nique s (e.g., iGuaita et^l2010l : iBerrv etani2012t IVargas et aP 
l2014l) which bracket strong emission lines at z ^ 2. All these ap¬ 
proaches select star-forming galaxies, and yield insight into the 
star-forming properties of these galaxies, but they have differing 
selection bia ses, and therefore, these samples do no t completely 
overlap (see iLv et al.l l20 1 it [Haberzettl et al J l20 1 21 for details). 
Therefore, it is essential to apply a well defined selection lead¬ 
ing to the identification of a population with a broad range of 
galaxy properties, which can then be compared to similar galax¬ 
ies at higher redshifts. 


The comparison of stellar populations of (strong) Lya emit¬ 
ting galaxies (LAEs), primarily selected based on the NB imag¬ 
ing technique, and non-LAEs (with weak or no Lya-emission), 
primarily selected based on photometric colors that mimic a 
break (i.e., Lyman break galaxies, LBGs), ha s yielded diverse 
and inconclusive results. iGawiser et al.l (l2006l) found that NB- 
selected LAEs at z-3.1 are less massive and less dusty com¬ 
pared to continuum selected LBGs suggesting that LAEs repre¬ 
sent an early stage of an evolutionary sequence where galaxies 
gradually beco me more massive an d dusty due to mergers and 
star formation (IGawiser et alJl20()^ . This conclusion is consis¬ 
tent with high specific star formation rates (SSER, SE R per unit 
mass) of NB-selected LAEs found bv iLai et51 ( 20081) relative to 
LBGs at the same redshift. On the other hand, Einkelstein et al.l 
(l2009l) found a range of dust extinctions in a sample of 14 NB- 
selected LAEs at z - 4.5, which is not consistent with LAEs be- 
ing the first g alaxie s in the evolutionary sequence. Similarly, 
iNilsson et alJ (1201 If) used 171 NB-selected LAEs at z-2.25 
and concluded that the stellar properties of LAEs are different 
from those at higher redshift and that they are diverse. They 
also believe that Lya selection could be tracing different galax¬ 
ies at_differeny;edshift^which_2S_consistent with the findings 
of lAcauaviva et all (l2012l) . iLai et al.l (l2008l) found that IRAC- 
detected LAEs are significantly older (age > 1 Gyr) and more 
massive (mass ~ 10'° M©) compared to IRAC-undetected LAEs 
at z-3.1, which lead them to suggest that the IRAC-detected 
LAEs may be a lower-mass extension of the LBG population. 
These studies show heterogeneity of NB-selected LAEs, and 
comparison between these objects and continuum selected LBGs 
continues to be the subject of many new studies. 


Various authors have compared stellar population prop¬ 
erties o£^^AEs_md_non-UAEs_Jn_(JY flux-limited samples 


(e.s.. Shanlevetal. 20011:lErb etalJ20061: Pentericci etal.1 

20071; 

iReddv etal.ll2008l:lKornei etalj|2010l). Shanlevetal.l(2001 

)gen- 


erated rest-frame UV composite spectra of LBGs at z - 3 divid¬ 
ing subsamples in ‘young’ (f < 35 Myr) and ‘old’ (f > 1 Gyr) 
galaxies. They found that younger galaxies have weaker Lya 
emission, are more dusty, and less massive c ompared to older 
galaxies. On the other hand, the analysis of lErb et"^ (l2006l) 
of the composite UV spectra for ~60 SEGs at z - 2 concluded 
that, on average, objects with lower stellar mass had stronger 
Lya emission line than more massive objects. Using a sample 
of 14 UV-selected galaxies at z - 2-3 with Lya equivalent width 
(EW) > 20A, iReddv et"^ (l2008l) found no significant difference 
in the stellar populations of strong Lyg-e mitters compared to the 
rest of the sample. At higher redshifts, iPentericci et al.l (l2007l) 
found that, in general, younger galaxies at z - 4 showed Lya in 
emission while older g alaxie s showed Lya in absorption. More 
recently. iKornei et alJ (1201 Ol) used 321 LBGs at z-3 to study 
the relationship between Lya-emission and stellar populations. 
Based on their analysis, they concluded that objects with strong 
Lya-emission are older, lower in SER, and less dusty compared 
to objects with weak or no (Lya in absorption) Lya-emission. 
The results of these studies emphasize that the exact relation be¬ 
tween stellar populations and Lya-emission is still not yet fully 
understood. 

In this paper we use spectroscopic data from the VIMOS 
Ultra-Deep Survey (VUDS), to study the spectral & photomet¬ 
ric properties of a large sample of faint SEGs at 2<z<2.5 
(Zmedian -2.3). VUDS recently obtained data for ~10000 galax¬ 
ies over an area of ~1 deg^ using the 8.2m Very Large Tele¬ 


scope (VLT). The VUD S observations were performed in the 
well-studied COSMO S (IScoville et al.l 2007 : ntoekemoer et al.l 
l2007l) and ECDES (iGiavalisco et al. 


20041 : iRix et all 120041) 
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fields now partly covered with HST/WF C3 by CANDELS 
dGrogin et alj 20Hl; Koekemoer et al. 201 Ih, and i n the VVDS- 
02h field ( Le Fevre et al.l 120041 : iMcCracken et'aDl2003h . These 
fields have extensive multi-wavelength data. We extract a sam¬ 
ple of 854 galaxies with confirmed spectroscopic redshifts at 
2<z<2.5, uniformly selected from their continuum properties 
with Iab ^ 25 mag. The goal of this paper is to study the phys¬ 
ical properties of these SFGs which include galaxies with and 
without Fya emission. 

This paper is organized as follows: In ^ we summarize the 
VUDS observations, and discuss 2<z<2.5 sample selection as 
well as its basic properties. In |3] we fit observed rest-frame UV 
VUDS spectra at 2 < z < 2.5 to measure the UV spectral slope, 
and discuss its correlation with the Es(B-V), stellar mass, and 
UV absolute magnitude. In ^ we discuss differences in stellar 
population properties between galaxies with and without Fya 
in emission. In ^ we investigate correlations between derived 
physical parameters (UV absolute magnitude, SFRs, Es(B-V), 
stellar mass, and UV spectral slope) as a function of rest-frame 
Fya EW for SFGs and their implications on our understanding 
of these galaxies. Results are discussed in |6] and we conclude 
with a summary in ^ 

Throughout this paper, we assume the standard cosmology 
with Q.„-03, Q.a-0.7 and Hq—IO km s ' Mpc '. This corre¬ 
sponds to a look-back time of ~10.4 Gyr at z -2. Magnitudes 
are given in the AB system (lOke & Gunnlll98^ . 


2. Data and Sample Selection 

The VUDS observations (iFe Fevre et al.ll2015l) were done using 
the low-resolution multi-slit mode of VIMOS on the VET A to¬ 
tal of 15 VIMOS pointings (~224 arcmin^ each, ~1 deg^ total) 
were observed with both the FRBFUE and ERRED grisms cov¬ 
ering the full wavelength range from 3650A to 9350A in three 
fields (ECDFS, VVDS-02h, COSMOS). The total exposure time 
per pointing was 14h in each grism. This paper is based on ~80% 
of the data, which is the amount processed at the time of writ¬ 
ing this paper. Details of these observations and the reduction 
process are described in iFe Fevre etal.1 (l2015l) . 

2.1. Photometry 

VUDS targeted three well-studied extragalactic fields with 
extensive mul t i-wave length photometry. Details are given in 
iFe Fevre et al.l (l2015h . Here, we briefly summarize the most rel¬ 
evant imaging observations for these three fields. 

The COSMOS field was observed wi th HST/ACS in 
the F814W filter (IScoville et al.l l2007t iKoekemoer et al.l 
l2007l) . Deep ground based imaging includes observations in 
u*,B,V,g^,r^,P,z^ bands from Subaru and CFHT. Details 
about these and other publicly available multi-wavelength 
observations are available at COSMOS website^]. The Ultra- 
Vista survey is acquiring very deep near-infrared imaging in 
the YJH K bands using the VIR CAM camera on the VISTA 
telescope (iMcCracken et al.l[2012l) . while deep Spit zer/IRAC ob¬ 
serva tions are available through the SCOS MOS dSanders etakl 
2QQ7b program. The C ANDEFS survey dGrogin et al.l 120111 
Koekemoer et al.l [201 ll) also provided WFC3 NIR photometry 
in the smaller, central part of the COSMOS field. 

The ECDFS field has been the target of many multi¬ 
wavelength surveys in the last decade. The central part of 

' http://cosmos.astro.caltech.edu/data/index.html and 

http://cesam.lam.fr/hstcosmos/ 


the fi eld, covering ~ 160 arcmin^, has HST/ACS observa¬ 
tions dGiavalisco et al.l |2004|) combined with the recent CAN- 
DEFS HST/WFC 3 observations in the near-infrared bands 
dKoekemoer et al.l 1201 ill . The ECDFS field is covered by deep 
UBVRIzJ HK imaging as part of MUSYC and other surveys in 
this field dCardamone et alj|2010l and references therein). The 
SERVS pro gram obtained mediu m-deep observations in 3.6//m 
and 4.5jum dMauduit et al.ll20l^ , which complement those ob¬ 
tained by the GOODS team (PI: M. Dickinson) at 3.6, 4.5, 5.6, 
and 8.0jum. 

The VVDS-02h field has dee p optical imaging from CFHT 
in the u*,g',r’,i',z' bands (e.g.. ICuillandre etakl I2OI2I) as part 
of the CFHT Fegacy Survey. Deep near-infrared imaging in 
the JHKs bands has been obtained with WI RCam on CFHT a s 
part of the WIRCam Deep Survey (WIRDS: lBielbv etal.ll2012h . 
The deep near-infrared imaging data (3.6jum and 4.5jum) has 
been obtained with Spitze r/IRAC as part of the S ERVS program 
re.g..lMauduit etaDl2012l) . We refer the reader to iFemaux et al.l 
( 20 Mall for a detailed description of multi-wavelength imaging 
data in VVDS-02h field. 


2.2. VUDS Target Selection 

The primary selection criterion for galaxies in the VUDS pro¬ 
gram was photometric redshifts (Zphot + ter >2.4), which are 
accurate to within 5% errors for these well-studied fields (e.g., 
Illbert et ^ 120091: ICoupon et akl 1200^ iDahlen et al.l l20ldl) . For 
high redshift (z > 2.5) targets, this photometric redshift selection 
was supplemented by color criteria (e.g., FBG/ugr, FBG/gri, 
FBG/riz), if the objects satisfied these color criteria but were 
not selected from the primary photometric redshift criterion. The 
fraction of targets selec ted for each criterion a t different redshifts 
is shown in Table 2 of iFe Fevre et al.l (l2015l) . At 2 < z < 2.5, all 
spectroscopic targets were selected solely based on the photo¬ 
metric redshift criterion. Therefore, the targets for the VUDS 
program include a representative sample of all star-forming 
galaxies at a particular redshift within a given magnitude limit 
Uab ^ 25 mag, with some galaxies as faint as iab ~ 2.1 mag). It 
is important to note that the target selection is not based on any 
particular emission line of a galaxy but on continuum magni¬ 
tude. A detailed discussion about the target selection, reliability 
of the redsh ift measurements and corresponding quality flags is 
presented in lUe Fevre etal.1 (l2015l) . 

2.3. Sample Selection 

For this paper, the primary selection criterion for the sample 
of SFGs from the VUDS spectroscopic data is the redshift. All 
objects between z-2 and z = 2.5 are selected in the final sam¬ 
ple, keeping only the best reliability flags (2,3,4,9) — which 
gives very high probability (75 -85%, 95-100%, 100%, 80%, re¬ 
spectively; see lUe Fevre et al.l uMsI for details) for these red¬ 
shifts to be correct. The total number of galaxies selected at 
2 < z < 2.5 for the current work is 854. Based on a simple nearest- 
neighbor matching to the most recent catalogs in each fi eld (e.g., 
IChiappetti etaH2005tlEivis etal.ll200^[Xue etal.l20llh . we find 
that our sample contains minimal contamination from X-ray 
AGN (<1%), which suggests that our sample is comprised al¬ 
most exclusively of galaxies without powerful X- ray AGN ac¬ 
tivity. We also checked for obscured AGN using IPonlev et al.l 
(I 2 OI 2 I) Spitzer/IRAC criteria, and found that ~3% of galaxies 
in our sample could host AGN based on this criteria. Some of 
these IRAC-selected AGN are also X-ray AGN, so after remov- 
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ing these overlapping AGN we end up with ~2% of additional 
AGN candidates using the Donley et al. IRAC criteria. This con¬ 
firms that our SFGs sample has minimum contamination (<3%) 
from AGN identified based on their X-ray emission and IRAC 
colors. This does not rule out the presence o f other types o f AGN 
or AGN-like activi ty in our sample (e.g., ICimatti et al.l 120131 : 
iLemaux et alJ2010h . but such phenomena do not, typically, dom¬ 
inate the UV/optical/NlR light of their host galaxies, which is the 
important point for our analysis. _ _ 

The Le PHA RE software package (lArnouts et alJ 1 19991 : 
lllbert et ^l2006h was used to fit the broad-band observed Spec¬ 
tral Energy Distributions (SED s) with synthetic stellar popu- 
lation models. We use BC03 dBruzual & Charlol l2003h tem¬ 
plates to generate a set of stellar populatio n models assum¬ 
ing a Chabrier initial mass function (IMF; IChabrieij l2003h . 
and fixed the spectroscopic redshift, while varying metallic- 
ity (Z = 0.4 and 1 Zq), age (0.05 G yr <t< tn), d ust extinction 
(0 < Es(B-V) < 0.7 mag) — using a lCalzetti et al.l (l2000h atten¬ 
uation law or, in some cases, Calzetti law -H 2175A bump — 
and e-folding timescale (t = 0.1-30 Gyr, i.e., from instantaneous 
burst to continuous star formation) for a star-formation history 
(SFFI)ocexp(-t/T). We also included two delayed SFH models 
with peaks at 1 and 3 Gyr. The Le PHARE code assumes the 
iMadaul (1 19951) prescription to estimate inter-galactic medium 
(IGM) opacity. From the best-fit model, we estimate stellar mass, 
dust extinction Es(B-V), SFRs, and SSFR for each galaxy. The 
UV absolu t e mag nitudes are measured using the prescription of 
lllbert et ^ (120051) and are not corrected for the internal dust ex¬ 
tinction. 

We emphasize that SFRs and Es(B-V)/dust inferred from fit¬ 
ting BC03 models to the SED of galaxies give only a rough es¬ 
timate of these physical parameters. To understand the effect of 
medium-band filters as well as other input quantities on the best- 
fit SED parameters, we compared our SER and Es(B-V) mea¬ 
surements with two different surveys that use medium band pho¬ 
tometry: (1) NEWEIRM Medium Band Survey (NMBS) in the 
COSMOS field dWhitaker etal.ll2012ah and ( 2) COSMOS SEP 
measurements using 30 photometric bands (lllbert et al.l l2009l) . 
which includes broad-band as well as intermediate and narrow 
bands in the optical and NIR. Eor NMBS data, we have a small 
number (~20) of galaxies matching both catalogs. Eor these 
galaxies, we find that, on average, SERs differ by ~0.5 dex (fac¬ 
tor of ~3), and A„ (E,(B-V)) differ by < 0.2 (< 0.07). Eor COS¬ 
MOS 30-band data which has larger sample (~300) of galaxies, 
we find that, on average, SERs differ by ~0.2 dex (factor of ~2), 
and A,, (Es(B-V)) differ by <0.15 (<0.05). It is crucial to note 
that NMBS uses slightly different input set-up for SED-fitting 
e.g., slightly different SFHs, different choice of extinction laws, 
different metallicity range, use of photometric redshifts, com¬ 
pared to what we use, so these choices of input parameters will 
also affect SERs and dust measurements. Therefore, the larger 
difference between NMBS and our SED measurements is not 
only due to photometric bands. As for COSMOS 30-band, we 
kept all the input parameters same as in the broad-band photome¬ 
try SED fitting but changed only the input photometric bands, so 
other effects do not play any role, and we see much smaller dif¬ 
ferences. These comparisons show that the uncertainties in SERs 
and dust measurements due to the choice of photometric bands 
and other input SED parameters could be as high as factor of ~3 
for SERs and ~0.2 mag for A„ (or ~0.07 for ETB -V)), which is 
typically expected for the SED-fitting process (e.g.. lMuzzin etalJ 
120101: iMostek et alj|2012|) . It is important to mention that these 
uncertainties are inherent in the SED fitting process and are not 
unique to this study. Moreoyer, this study is based on the rela- 


tiye measurements of three different galaxy populations and as 
long as there are no drastic differences in the IME or SEHs of 
these populations, these systematic uncertainties does not play a 
significant role in the conclusions of this paper. 

The contribution of major emission lines in different filters 
has been included in the models of the Le PHARE code. We 
use these templates with emission lines to carry out SED fit¬ 
ting. Neglecting emission lines during the SED fitting process 
can oyerestimate t he best-fit stellar ages and stellar masses by up 
to ~ 0 .3 dex (e.g . . ISchaerer & de BaiTosll2009t iFinkelstein et ai.l 
1201 ll : lAtek et akl 1201 ih . The Le PHARE code accounts for the 
contribution of e mission lines with a simple recipe based on the 
iKennicutl (Il998h relations between the SER and UV luminosity. 
Hot and [Oil] lines. The code includes the Lya, Ha, Hff, [Oil], 
OIII[4959] and OIII[5007] lines with different line ratio s with 
respect to the [Oil] line, as described in lllbert et al.l (l2009l) . 

2.4. Sample Properties 

Eigure [1] shows distributions of spectroscopic redshifts and far- 
UV (1500A) absolute magnitudes (M 1500 ) for the full SEG sam¬ 
ple at 2 < z < 2.5. The median M 1500 of the sample is -20.22 mag, 
while the median redshift is 2.320. The figure also shows the 
eyolying characteristic magnitude M* a s a function of the red¬ 
shift, as deriyed from the compilation of lHathi et al.l (l2010l) . The 
M*+l mag range coyers most of our data which implies that we 
are probing luminosities around L* at these redshifts. The shape 
of the redshift distribution is dictated, in part, by the selection 
function used to select target galaxies for VUDS. The main se¬ 
lection criteria under which primary targets are selected, is Zphot 
-H Icr > 2.4. Therefore, at lower redshifts (2 < z < 2.4) we are on 
the tail of the probability distribution of photometric redshifts, 
and this explains the decreasing number of galaxies at z < 2.4. 
This bias equally affects all galaxies — galaxies with and with¬ 
out Lya in emission — so the general results in this paper, mak¬ 
ing relatiye comparisons of the two populations, are not affected 
by this selection bias. The second reason for the small number 
of galaxies close to the lower limit of our redshift range could be 
the sensitiyity of the VIMOS instrument which starts to decrease 
dramatically below 3800A. This means that we might not be able 
to properly identify continuum breaks, important features to se¬ 
cure the redshift, below redshift z < 2.15. In this redshift range, 
we might be biased in fayor of strong Lya emitters as we could 
still identify those eyen at z<2.15. To check this sample selec¬ 
tion effect, we cut our sample at redshift z = 2.15. By doing so, 
we remoye ~10% of galaxies from the full sample, which is yery 
small compared to the large sample size of VUDS galaxies. This 
small change in number of galaxies does not change our final 
results. Therefore, the results presented in this paper are not af¬ 
fected by this selection bias. 

Eigure |2] shows distributions of stellar mass and SER for the 
full SEG sample at 2 <z<2.5. The median stellar mass of the 
sample is ~4xl0^ Mq and the median SER is ~20 Mp per year. 
We sh ow the ‘main sequence’ of SEGs as defined by lDaddi et al.l 
(I 2 OO 7 I) at z - 2, and the best-fit relation for our sample with 
the corresponding scatter (~0.4 dex). We keep the logarithmic 
slope fixed to the yalue from Daddi et al. relation because the 
redshift probed is yery similar for both studies. If we fit the ob- 
seryed SER-stellar mass relation without keeping the slope fixed, 
then we get the logarithmic slope of ~0.66 which is flatter than 
the slope of 0.90 obtained by Daddi et al. at similar redshifts. 
Eigure |2] shows that our sample has a large scatter in the stel¬ 
lar mass-SER plane coyering a large range in SER and oyer two 
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Fig. 1. Redshift versus UV absolute magnitude (M 1500 ) for spectroscopically confirmed SFGs at 2<z<2.5 from VUDS. The median M 1500 of 
the sample is -20.22 mag, the median redshift is 2.320. The quoted uncertainties are the errors on the median (1.253 x cr/ -y/N^) estimated 
from the Icr dispersion of the distribution and the total number of galaxies, where ± Itr corresponds to the range between the 16th and the 84th 
percentile values of the distribution. The red solid line indicates the evolving M* (characteristic magnitude) as a function of redshift as derived 
from the compilation of lHathi et al.l ll2010h . The dashed red lines indicate M*±l mag. The black solid curve on the histogram is the kernel density 
estimation (KDE) of the distribution, which is a smoother non-parametric density estimator compared to the histogram and is not affected by the 
bin-size or the end points of the bins. The density of points is color-coded as shown in the color-bar. 


orders of magnitude in stellar mass. It is important to note that 
our best-fit relation (keeping the logarithmic slope fixed at 0.90) 
between SFR and stellar mass is ~0.2 dex off from the Daddi’s 
SF main sequence. The Daddi et al. sample is at slightly lower 
redshift compared to our sample, so we also show in F i gure |3 
the best-fit main sequence relation from IWhitaker et al.l (12014h 
and ISchreiber et al.l (l2015h . The redshift probed by these stud¬ 
ies is very similar to ours. We still find a similar offset (~0.2 
dex) between our main sequence relation and the relation from 
these studies, which predicts slightly lower normalization i.e., 
lower SFR at a given stellar mass. This difference could be due 
to the fact that we estimate SFRs based on the SED fitting pro¬ 
cess which has wide variety of SFHs and it gives slightly differ¬ 
ent SFR compared to SFR(UV) or SFR(UV-i-IR). When we esti¬ 
mate SFR using UV luminosity (corrected for dust using the UV 
slope), we estimate average SFR to be ~0.2 dex lower compared 
to SFR(SED). Therefore, this could be a reason for ~Q.2 dex 
offset we see in the main sequence relation compared to other 
studies, as they use SER(UV) or SER(UV-i-IR) in their relations. 
By selection, and confirmed by NUV - r - J plot, our galaxies 
are not very dusty and hence, not many galaxies have IR/EIR de¬ 
tections. In future papers, we will investigate IR SERs for these 
galaxies based on their 24//m and Herschel observations. 

The scatter (total scatter obtained from the fitting process) 
in the SER- stellar mass correl ation is ~QA dex, which is higher 
than that of iDaddi et al.l (l2007h by ~Q.2 dex. We believe that this 
large scatter is caused by the choice of wide range of SEHs in the 


SED fitting process . Various studies (e.g.. iMcLure et al1l201 ll: 
ISchaerer et al.ll201^ have shown that a broader range in SEHs 
in the SED fitting process could have strong implications on the 
scatter of the main sequence relation. Considering that our sam¬ 
ple has similar redshift (z - 2) as Daddi et al. and other studies, 
these differences in the main sequence relation emphasize that 
the sample properties as well as measurement techniques play a 
significant role in the determination of the co rrelation between 
SER and stellar masses of galaxies (see, e.g., IPforretakllMl 
for uncertainties in SED fitting process). To underscore the im¬ 
portances of these differences, we point out few subtleties be¬ 
tween the Daddi et al. relation and ours. Our sample was se- 
lected based on th e photometric redshift technique, while the 
iDaddi et^ (l2007h sample was selected based on a single BzK 
color technique, highlighting that sample selection likely plays 
a significant role in determining the parameters of the stellar 
mass-SER relation. We use different star-formation histories in 
our S ED-fitting technique that estimates SER and stellar masses, 
while IDaddi et^ (l2007l) obtain SERs directly from UV lumi¬ 
nosity and their stellar m asses are based on a single color tech¬ 
nique (IDaddi et al.ll2004l) . This difference in method also likely 
influences the measurements of these stellar parameters and thus 
the shape & normalization of the stellar mass-SER relation. The 
detailed discussion about various aspects of this relation is be¬ 
yond the scope of this paper, therefore, for in-depth study on this 
relation and its e volution with redshift, we refer the reader to 
iTasca et ^ (l2015l) . and for more details on the effect of different 
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Stellar Mass [Log ( M^)] Number 

Fig. 2. SFR versus stellar mass for SFGs at 2<z<2.5 in VUDS. The full sample has median stellar mass of 10®'®'* ~4xl0® Mq and median 
SFR of 10*'^^ ~20 Mq per year. The quoted uncertainties are the errors on the median (1.253 x cr/ estimated from the Itr dispersion of 

the distribution and the total number of galaxies, where ±lcr corresponds to the range between the 16th and the 84th percentile v alues of the 
distrib ution. The black solid curve is the KDE of the distribution. The black solid line shows the main sequence of SFGs as defined bv iDaddi et al.l 
ll2007h at z ~ 2, and the solid red line is the best-fit relation (with the slope fixed to the value from Da ddi et al.) for our sam ple with the total 
scatter (~0.4 dex) shown by red d ashed lines. The black dashed line is the main sequence relation from IWhitaker et ^ (l2014ll . while the black 
dot-dashed line is the relation from ISchreiber et £0 ( 120151) . The blue dot-dashed lines show limits on SSFR for these galaxies. The density of points 
is color-coded as shown in the color-bar. 


SFHs on the SFR-stellar mass correlation, we refer the reader to 
Cassara et al. (in prep). Both these investigations are based on 
the same VUDS dataset, but over a larger redshift range. 

3. UV Spectral Slope {J3) 

The UV spectral slope (J3) is determined from a p ower-law fit to 
the UV continuum spectrum (ICalzetti etal.ll994l) . /r cc where 
is the flux density per unit wavelength (ergs s“' cm“^ A“'). We 
use the VUDS spectra of SFGs at 2 < z < 2.5 to estimate their UV 
slope by fitting a linear relation to the full wavelength range be¬ 
tween rest-frame 1400 and 2400A in the observed (Log(/i) ver¬ 
sus Log( /;)) spectrum. The rest-frame wavelength range (1400- 
24 00A) includes 7 out of 10 spectral fitting windows identified 
by ICalzetti et al.l (Il994l) to estimate the UV spectral slope for lo¬ 
cal starburst galaxies, and takes full advantage of the spectral 
coverage of the VUDS spectra at these redshifts. Figure |3] shows 
four examples of UV spectral slope fitting on the VUDS spectra. 
Before fitting, the spectra were cleaned by masking out known 
lines and removing any spurious noise spikes. The fitting was 
done using the IDL routine UNFIT, which fits the data to the lin¬ 
ear model by minimizing the;^'^ statistic. The running dispersion 
values measured on the clean spectra are used as measured er¬ 
rors, and the routine uses these errors to compute the covariance 
matrix. In Figure[3l the solid red line is the best-fit spectroscopic 


slope while dashed red lines show Icr uncertainties associated 
with the best-fit relation. The formal statistical uncertainty in the 
best-fit ySspec. which is measured over a large number of points, 
is smaller than the uncertainty in photometric /? measured us¬ 
ing few photometric points. For comparison, we also measure 
the UV spectral slope directly from the broad-band photometry 
which encompasses a similar rest-frame wavelength range. We 
used BVRI bands for the ECDFS held, g^Vr^i^ for the COS¬ 
MOS held, and g'r'i'z' bands for the VVDS held. These photo¬ 
metric bands cover a rest-frame wavelength range from ~1200A 
to ~3000A, depending on the redshift. Figure |3] also shows the 
photometric magnitudes of the broad-band Alters which are used 
to measure the photometric UV slope, and the corresponding 
best-fit photometric UV slope. Figure |4] shows the UV spectral 
slope (yS) distribution as measured from the VUDS spectra (ySspec) 
and the broad-band photometry (ySphot)- The median ySspec is - 
1.36 + 0.02, and the median/3phot is -1.32 + 0.02. The quoted un¬ 
certainties are the errors on the median (1.253 x cr/ It is 

worth noting that even if we use two or three photometric bands, 
covering the similar wavelength range as four bands, we get sim¬ 
ilar median values for the photometric UV slope of the whole 
sample, which is in the range of-1.40+0.10. We also measured 
UV slopes by convolving the spectra with the photometric Alter 
curves, and found that the median value is -1.31 +0.03, which 
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Fig. 3. Examples of UV spectral slope (/?) fitting on the VUDS spectra (shown in grey). The spectra between rest-frame 1400A to 2400A (between 
dashed vertical lines) are fitted with a linear relation (solid red lines) between Log(/l) vs Log(Fa) after masking out known lines and removing 
spurious noise spikes. The dashed red lines are Icr deviations in the best-fit relation. The blue filled circles are the photometric magnitudes from 
the broad-band filters (grey curves) used to measure the photometric j8, while the blue dashed line is a linear fit between these photometric data 
points to estimate the photometric UV slope. 


is very similar to the median photometric and spectroscopic UV 
slopes. 

In Figure |4], and in subsequent figures of this paper, we will 
calculate statistical significance for the distributions and correla¬ 
tions in the following two ways: (a) for histograms, we use KS- 
statistics on the distributions. The probability of KS-statistics 
(Pks) measures whether two sets of data are drawn from the 
same parent distribution or not. Small values (close to zero) of 
Pks show that the cumulative distribution functions of the two 
sets of data are significantly different. This is measured using 
the IDL routine KSTWO, and (b) for correlations, we use the 
Spearman correlation test. Values of the Spearman correlation 
coefficient (p) close to +1 imply strong monotonic correlation, 
while values close to 0 imply no significant correlation between 
two variables. The significance level is denoted by Psc, where 
small values of Psc imply significant correlation. It is important 
to note that for a large sample size, very small differences (or 
small p values) in two sets of data will be detected as highly 
significant. When a statistic is significant, it simply means that 
the difference is real. It does not mean the difference is large or 
important. The Spearman correlation coefficient is measured us¬ 
ing the IDL routine R_CORRELATE, and is measured from all 
galaxies as well as from the median-binned values.This routine 
can also give the value of the sum-squared difference of ranks, 
and the number of standard deviations (cr) by which it deviates 
from its null-hypothesis expected value. 


In Figure |4] (left panel), the median values of /3spec and ySphot 
are very similar and differ by only ~0.04. This difference is 
not significant considering that the errors on the median val¬ 
ues is +0.02. However, the full distributions of ySspec and ySphot 
are marginally different as confirmed by Pks ~ 0.00057, which 
implies that a null hypothesis (similar distributions) is rejected 
at ~3cr level. To investigate the difference between individual 
ySspec andySphot. we look at four examples shown in Figure[3] Fig¬ 
ure [3a) shows an extreme example in which the difference be¬ 
tween two J0S is ~0.5, thoughySphot is within the Icr uncertainties 
of ySspec- Figure 13b) and Figure |3c) show examples where ySphot 
is just outside the Icr uncertainties in ySspec and we find that both 
these p values differ by ~0.2-0.3. Figure|3d) shows an example 
where the difference between the two jSs is very small (~0.1) and 
well within the measurement uncertainties. 

The middle panel of Figure |3 shows the direct compari¬ 
son between individual ySspec and /?phot, while the bottom plot in 
the middle panel shows the difference between ySspec and ySphot 
as a function of ySspec- The best orthogonal fit relation between 
ySspec and ySphot is flatter than 1-to-l relation, and is given by 
ySphot = -0.43 + 0.63 x ySspec- The total scatter in the jSspec vs 
ySphot relation could be as high as ~0.5. The uncertainties in both 
P measurements could be the main cause of the flatter slope (< 1) 
and the large scatter. The average measurement error in ySspec is 
~0.08, while it is ~0.15 for ySphot- Depending on the object and 
its redshift, it is possible to have additional factors affecting p 
measurements of individual galaxies, as discussed below. 
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For ySspec, there are two potential sources of uncertainties. 
First, the bluest part of the spectra is sensitive to various at- 
mospheric and galactic corrections mad e to the spectra (see 
iLe Fevre et al.ll20r^ iThomas et alJl20l4 for details), and even 
though these corrections are carefully applied to minimize the 
flux loss, there are uncertainties in the bluest part of the spec¬ 
tra. The UV slope measurement is sensitive to the changes in 
the anchor blue region of the spectra and hence, this could affect 
the UV slope measurements for some objects. Second, as seen 
in Figure [3a), the spectroscopic fi is fitted in such a way that it 
avoids the more noisy red part of the spectrum to minimize its ef¬ 
fect on the best-fit slope, which could make /Ispec steeper or bluer 
compared to jSphot- While this fitting effect is kept at the mini¬ 
mum level, any small difference in the fitting of the red part of 
the spectrum could lead to slightly bluer ySspec compared to jSphot- 
These uncertainties in ySspec are estimated to be small (<0.1) but 
could vary from object to object. 

ForySphot, the wavelength range over wh ich /? is measu r ed has 
a great influence on the resulting values. ICalzetti etaP (Il994ll 
defined a wavelength baseline (1250A to 2600A) to measure 
j3 for galaxies at lower redshifts and this r ange is widely used 
for all measurements. But as discussed in ICalzettil (l200lh . the 
UV slopes fitted with the longer wavelength baseline are always 
bluer (more neg ative) compared to the UV slopes fr om shorter 
baseline (see also lMeurer et alJl999HTalia etalJ20T3) . The main 
reason is that the slopes fitted with the longer wavelength base¬ 
line cover the 2300-2800A range which has a large decrement 
due to a large number of closely spaced Fell absorption lines and 
this low continuum gives them steeper/bluer slopes. Our jSphot 
is measured using broad photometric bands, and depending on 
the exact redshift of a galaxy the wavelength range covered by 
these bands could be as large as ~1200A to ~3000A. This range 
includes the wavelength range (1400A to 2400A) used to mea¬ 
sure jSspec, but ySphot has a longer wavelength baseline for many 
galaxies which could lead to bluer jSphot for these galaxies. The 
other factor affecting the jSphot measurement is the contribution 
of strong Lya emission to the flux of the bluest band used for 
measuring jSphot- For a small number of strong Lya emitters, it is 
likely that the bluest band used in the jSphot measurement is af¬ 
fected by the strong Lya emission line. To check this effect, we 
measured j8s for galaxies with no Lya emission and for galaxies 
with strong Lya emission. We find that the difference between 
the median values of /Iphot and jSspec is smaller (ySspec bluer by 
~0.05) for galaxies with no Lya emission compared to galax¬ 
ies with strong Lya emission, where on average ySphot is bluer 
compared to ySspec by ~0.13. Hence, strong Lya line in the bluest 
photometric band could lead to higher flux in that band which in 
turn could give steeper/bluer jSphot- This is confirmed by our esti¬ 
mation of Lya contamination to the bluest photometric band by 
comparing flux densities with and without Lya line in the broad¬ 
band. We found that the contamination in the broadband for the 
strongest Lya emission lines could be ~0.05-0.1 mag. This dif¬ 
ference in the magnitude of bluest band could cause ~0.1-0.2 
difference in the photometrically measured UV slope. These un¬ 
certainties in ySphot are estimated to be <0.2, but could vary from 
object to object. 

In summary, the average jSspec and jSphot are very similar for 
the VUDS SFG sample, and the uncertainties discussed above 
could affect individual measurements, as shown in the middle 
panel of Figure |4l The detailed comparison between individual 
jSs, and quantifying various sources of uncertainties for such a 
large sample, requires an in-depth study, which is beyond the 
scope of this paper. 


The median ySspec value of the whole SFG sample, - 
1.36+0.02 with Icr scatter of ~0.5, is consistent with the UV 
slope estimate s obtained at sim ilar redshifts using ground- 
based spectra. iTalia et all (1201 2ll obtained ySspec of -1.11+0.44 
(rms) for 74 SFGs a t z^2 using VLT/FORS2 spectra, while 
iNoll & Pierinil (l2005h used FORS spectroscopic data of 34 UV- 
luminous galaxies at 2<z<2.5 and found /Ispec for the sub¬ 
sample of galaxies to be -1.01+0.11. Both these fi values are 
roughly consistent with our median /Ispec when considering the 
additional uncertainty introduced by the fact that both stud¬ 
ies used slightly different wavelength range (~1200-2600A, 
and ~1200-1800A, respectively), from each other & our own, 
to fit the UV slope. In addition, the median photometric UV 
slope obtained here, -1.32+0.02 with Icr scatter of ~0.5, is 
very similar to the average UV slopes measured at these red¬ 
shifts (2_<z_<2A}from HST photometry (e.g., -1.58 +0.27 (Icr) 
from Bouwens et alJl206^ and -1.71+0.34 (Icr) from iHathi et al.l 
12013 1. These p values — spectroscopic and photometric — are 
also in good agreement with the evolution of p with redshift 
which shows that the average at 2 < z < 2.5 is redder compared 
to P at higher redshifts (e.g Hathi et aDl2008bt iBouwens et al.l 
l2(mlFinkelstein etalJl20ia . 


3.1. p versus E^(B-V) 


The UV radiation is absorbed by dust and heated dust parti¬ 
cles re-emit absorbed energy in the mid- and far-infrared. The 
ratio of LiR-to-Luv is the direct measure of dust attenuation 
in a galaxy. A higher Lir/Luv ratio implies a larger dust con¬ 
tent. There is an excellent corre lation between Lir/Luv and the 
UV spectral slope at z ~ 2 (e.g., ISeibert et al.ll2002l iReddv et^ 
l2012h. which is in good agreem ent with the local measurements 
(e.g., iMeureretaklll^ fl^ . in that, an increase in Lir/Luv 
corresponds to redder UV slopes. For SFGs at 2 < z < 2.5 in the 
VUDS sample, the correlation between the UV slope and SED 
derived Es(B-V) is shown in Eigure |4| (right panel). A positive 
correlation is observed between p and Es(B-V) implying red¬ 
der UV slope for dustier (higher Es(B-V)) galaxies. We show 
the correlation of Es(B-V) with both, photometric and spectro¬ 
scopic, ySs. The statistical significance of this correlation is mea¬ 
sured by the Spearman correlation test. The Spearman correla¬ 
tion coefficient (p) is 0.25 (0.95) for all (binned) ySspec values and 
0.65 (1.00) for all (binned) ySphot values, which implies high sig¬ 
nificance (Psc ~ 0) for a strong positive correlation. The signifi¬ 
cance of p as measured by the number of standard deviations by 
which the sum-squared difference of ranks deviates from its null- 
hypothesis expected value is ~lcr for all ySspec values and ~18cr 
for all ySphot values. The correlation between ySphot and Es(B-V) 
is stronger, in part, because they are both estimated from the 
same photometry. The ySspec correlation is flatter compared to the 
ySphot correlation, which means that galaxies with bluer p have 
lower Es(B-V) for ySphot, and galaxies with redder p have higher 
Es(B-V) for ySphot, compared to jSspec- Also, it should be noted 
that within the large Icr scatter in the Es(B-V) values, the dif¬ 
ference between these two jSs is not large in all UV slope bins. 
The corresponding Aigoo (the dust absorption at 1600A) in Eig¬ 
ure |4| (right panel) is estim ated from Es(B-V) using the Calzetti 
law (Aieoo ~ lOxETB-VL iCalzetti eta^l.ll2000ll. We have added 
best-fit relations from three differ ent st udies ( lal^_et^ I^15t 
lOteo etal1l20l4 iBuat et al1l2012h . The Oteo et airi ^l4l) rela¬ 
tion is for galaxies at z ~ 2, w hile the Buat et al. (12012 ) relation 
is for galaxies at 1 < z < 2. The lTalia etal. ( 20151) relation covers 
a larger redshift range, 1 < z < 3. We do not see any trend (with 
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redshift) between these relations and ours. Our relation is flatter 
compared to their relations and the difference is largest at redder 
UV slopes. This difference between our relation and theirs could 
be due to the sample of SFGs used in these studies. All these 
studies use UV+FIR selected galaxies which are primarily dusty 
galaxies with FIR detection in 24//m and/or Herschel observa¬ 
tions, while we do only UV selection and do not require FIR 
detection for our sample. This means that our sample does not in¬ 
clude most of the dusty galaxies at these redshifts, which would 
make this relation much steeper wit h larger E(B-V) values and 
redder UV slopes. ITalia et al.l (l2015h finds that the differences in 
various samples could also lead to a large dispersion (~0.7) in 
the dust-UV slope relation. Our future studies will specifically 
look at IR/FIR detections for these galaxies to investigate their 
dust properties. 

In the remaining sections of this paper we use the spectro¬ 
scopic UV slope and refer to ySspec as fi, unless stated otherwise. 
We are using ySspec in this paper because they have smaller mea¬ 
surement uncertainties, and the trends/relations explored in this 
paper are very similar for both spectroscopic and photometric 
UV slopes. 


3.2. p versus Misoo 


A color-magnitude relation has previously been observed at var¬ 
ious redshifts. The interpretation of this relation depends on 
which color is used. If the col or is between the rest -frame UV 
and the rest-frame optical (e.g.. lPaDovich et^l2004t) and spans 
the 4000A/Balmer break, then it is strongly dependent on the 
stellar population age, while if the color is blue-ward of the 
4000A break (e.g., iFinkelstein etal.ll 2012 l) then it is more sen¬ 
sitive to the dust extinction than age. Here we use rest-frame UV 
colors (i.e., UV spectral slope, yS) to infer the color-magnitude re¬ 
lation for SFGs at 2<z<2.5. Figure |5] shows the spectroscopic 
UV slope versus UV absolute magnitude (Misoo). We estimate 
the slope, dpidMisoo, by fitting a linear model to the median- 
binned values, p - a + dpfdMisoo x M 1500 . We find the best¬ 
fitting solution to be dpidMisoo - 0.00+0.04. Essentially no cor¬ 
relation of p with absolute magnitude is observed, which could 
imply that p does not change with UV luminosity for ~L* galax¬ 
ies. The Spearman correlation coefficient (p) is 0.18 (-0.4) for 
all (binned) Misoo-jS values, which implies no correlation at a 
high significance (Psc = 0.60). The p for all galaxies shows high 
significance (Psc < 10 “^) for a strong correlation but it is mis¬ 
leading as it is strongly affected by a small number of galaxies 
with significant redder UV slope. 


Various studies have investigated the color-magnitude rela- 
tion at high redshifts (z > 2 ), and the results are inconclusive. 
iReddv et~m (l2008h found weak-to-no correlation between dust 
or UV color (related as shown in Eigure |4li and R mag for UV- 
selected galaxies at 1.5<z<2.6. Similarly, iHeinis etd] (l2013h 
through their investigation of far-infrared/dust properties of UV 
selected galaxies at 1.5 found that the average UV slope is 
mostly independent o f the UV luminosity. On the other hand, 
iBouwens et al.l (l2009h studied the relation between the UV slope 
and M[/v for a sample of t/-band dropout galaxies at z ^ 2.5 
and found a positive color-magnitude correlation. At higher red¬ 
shifts (z > 3), similar disagreements have been observed for sam¬ 
ples based on the Lyman break color selection, in the sense 
that, some studies find no (or very weak) trends between P 
and UV magn itude (e.g., Ono_e^ .l 2010 HFinkelstein etal.ll 20 f^ 
[Dunlop et al.l |20 1 2| : ICastellano et al.f 2012 ). while other studies 
find a strong color-magnitude relation (e.g., I Wilkins etal.ll 20 lil 


Bouwens,et^ 1^1^ iRogers etafll2014t iBouwens etal.ll2()T^ . 
Bouwens et al. (1201 4t) also argue that this relation could be de¬ 
scribed by a double power-law with a different slope for fainter 
magnitudes. The differences in these various studies could be 
due to many systematic and/or physical reasons, such as selec¬ 
tion of galaxy samples (especially at lower redshifts), UV magni¬ 
tude and/or p measurements, intrinsic scatter in p, dynamic range 
in UV luminosity. Our study relies on a UV selected sample 
with spectroscopic redshifts but it lacks a significant dynamic 
range in UV luminosity as seen in Eigure |5] A detailed study 
of a uniformly selected spectroscopic sample covering a large 
range in UV luminosity is needed to understand the true nature 
of the color-magnitude relation at high redshifts. At least for ~L* 
(~0.5L* - 3.0L*) galaxies in a uniformly selected VUDS sample 
at 2 < z < 2.5 there appears to be no correlation between p and 

Mi500. 

3.3. p versus Stellar Mass 

The relation between p and M 1500 could be affected by the bi¬ 
ases and/or differe nces in the way different authors measure and 
define M 1500 (e.g.. lEinkelsfein efal.1120121) . On the other hand, 
there is a well known correlation between M 1500 and the stel¬ 
lar m ass such that, stel l ar mass increa s es as brightness in creases 
(e.g.. lStarketalJl200^ISawi(^l2012l : [Hathi et alj|2013h . which 
we can exploit to better understand the dependence of p on 
Mi 5 oo. Therefore, we can compare p to the stellar mass as an 
alternate way to explore the dependence of p on Misoo- Eig¬ 
ure |5] shows the correlation between the spectroscopic UV slope 
and the stellar mass (M,). We estimate the slope, dpidLog(M,), 
by fitting a linear model to the median-binned values, p = 
a + dpidLog(MP) x Log(MP). We find the best-fitting solu¬ 
tion to be dpidLog{Mt) = 0.18+0.04. A stronger correlation of 
P with stellar mass is observed compared to the/l-Misoo correla¬ 
tion, which suggests a redder UV slope (more dust) for massive 
star-forming galaxies. This is not necessarily surprising as star- 
formation generally increases with stellar mass in star-forming 
samples and in turn, the dust content generally in creases with 
increasing star-formation (e.g.. lLemaux et al.l2014bh . The Spear¬ 
man correlation coefficient (p) is 0.11 (1.0) for all (binned) M*-y8 
values, which implies high significance (Psc ~ 10“^) for a strong 
positive correlation. A similar trend between P and ste llar mass 
is also observed at z > 4.0 (e.g.. lEinkelstein etal.ll2012h . 

The VUDS survey targets star-forming galaxies brighter 
than Iab ~ 25 mag, so there are survey limitation/incompleteness 
given the limited area of the survey, in the sense that, we do not 
target many high mass galaxies (Log(M,) > 10.5), and we miss 
many fainter star-forming galaxies with lower stellar masses 
(Log(M,) < 9.0) because of the magnitude limit. These survey 
limitations affect the range in stellar mass range probed in our 
study, and hence, could impact the exact nature of our /3-stellar 
mass correlation. An extensive study covering a large mass range 
is required to fully understand this M*-/? correlation. 

The Mi 5 oo-yS and M,-/? correlations are similar if we use 
ySphot instead of ySspec- For /3phot, the linear slope values for the 
Mi 5 oo-yS and M,-yS relations are 0.01+0.05 and 0.32+0.05, re¬ 
spectively. 

4. SFGs with and without Lyor Emission 

The galaxies selected by their (strong) Lya-emission, are becom¬ 
ing an important probe of galaxy formation, cosmic reionization, 
and cosmology. Lya emission is an important diagnostic of phys¬ 
ical processes in star forming galaxies, in particular at cosmo- 
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Fig. 4. (Left) UV spectral slope (/3) distribution as measured from the VUDS spectra (red histogram) and the broad-band photometry (blue 
histogram). The solid curve is the KDE of the distribution. The median fi values as well as KS-test probability (Pks) show very small difference 
in these /is. (Middle) /3spec versus ySphot- The solid red line is the best orthogonal fit relation between these two /is, while the dashed lines show the 
uncertainties in the best-fit relation. This figure also shows the difference (/SspeciSphot) as a function of/3spec. The dotted line is for = 0- The 
density of points is color-coded as in previous figures. Average uncertainties in these measurements are plotted in the hottom-right corner. (Right) 
The correlation between the UV slope and the SED derived Es(B-V) is such that the redder UV slope implies dustier [higher Es(B-V)] galaxies. 
The error bars in x illustrate the sizes of the bins, while the errors in y are ±lcr scatter (dashed error hars) corresponding to the range between the 
16th and the 84th percentile values within each bin, while the smaller solid error bars are the errors on the median values (1.253 x cr/ -y/N^). The 
red circles are for the spectrosc opic p, while the blue circles are for the photometric p. The black sm all-dashed line is t he dust-UV slope relation at 
1 ;£ z 5 2 from iBuat et al.l ll2012h . the black dot-dashe d line is the dust-U V slope relation at z ~ 2 from lOteo etalJ i2014ll . and the black long-dashed 
line is the dust-UV slope relation at 1 < z < 3 from iTalia et akl ( 120151) . The background grey points are the values for all galaxies (crosses are 
jSspec and open circles are ySphot)- The statistical significance from the Spearman correlation coefficient is shown to con firm the correlation . The 
corresponding Aisoo (shown on the right y-axis) is estimated from Es(B-V) using the Calzetti law (Aksoo ~ 10xE. (B-V): ICalzetti etalj20QOh . 



-18 -19 -20 -21 -22 

Mi 500 [mag] 


9.0 9.5 10.0 10.5 11.0 

Stellar Mass [Log ( Mg)] 


Fig. 5. Spectroscopic UV slope versus the UV absolute magnitude (left) and the stellar mass (right). The density of points is color-coded as in 
previous figures, while the red open circles are the median-binned values for these parameters. The error bars in r illustrate the sizes of the bins, 
while the errors in y are ±lcr scatter (dashed error bars) corresponding to the range between the 16th and the 84th percentile values within each 
bin, while smaller solid error bars are the errors on the median values (1.253 x cr/ The dashed red lines show the best-fit linear relations. 

We find no correlation (slope = 0.00±0.04) of p with the UV absolute magnitude, while we see a positive correlation (slope = 0.18±0.04) of p 
with the stellar mass (redder UV slope for higher mass galaxies). The statistical significance from the Spearman correlation coefficient is shown 
to confirm these results. 


logical distances, since Lya becomes the strongest emission line 
in the UV-optical window at redshifts z > 2.0. These galaxies — 
because of their observed physical properties at higher redshifts 
— are thought to b e the progenitors of p resent-day Milky Way 
type galaxies (e.g.. iGawiser et al.ll2()07h . though dissimilar re¬ 
sults from various Lya studies at different redshifts amplify the 
importance to better understand these galaxies. The SFG popula¬ 
tion in our sample is divided into Lya emitters and non-emitters 
based on their rest-frame Lya EW. 


The Lyg EWs for V UDS SEGs were measured as described 
in ICassata et al.l (l2015h . To summarize, we measured EWs of 
the Lya line manually using the IRAE splot tool, while the 
uncertaintie s in EW were estimated using the formalism from 
iTresse el'iD (Il999h . The median EW uncertainty varies as a 
function of the EW from +5A (for weak absorbers and emitters) 
to ±25A (for strong absorbers and emitters). We first put each 
galaxy spectrum in its rest-frame based on the spectroscopic red- 
shift. Then, two continuum points bracketing Lya are manually 
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marked and the rest-frame EW is measured for emission as well 
as the absorption lines. The line is not fitted with a Gaussian, but 
the flux in the line is determined by simply summing the pixels 
in the area covered by the line and the continuum. This method 
allows the measurement of lines with asymmetric shapes (i.e., 
with deviations from Gaussian profiles), which is the shape ob¬ 
served for most Lya lines. The interactive method also allows 
us to control the level of the continuum, taking into account de¬ 
fects that may be present around the line. This approach pro¬ 
duces more reliable and accurate measurements compared to an 
automated ‘impartial’ measurement because the visual inspec¬ 
tion is capable of accounting for the complex nature of the Lycr 
emission/absorption line. 

Figure |6] shows the rest-frame EW distribution for the Lya 
line. Here, positive EWs are emission lines, while negative EWs 
are for absorption lines. We divide the SFG population into 
three sub-groups based on their Lya EW. The galaxies that 
show no Lya in emission (EW < OA) are defined as SFGn, 
while the galaxies with Lya in emission, irrespective of its 
strength (EW > OA), are defined as SFGl- The third group is 
for strong Lya emitters (EW > 20A) called LAEs. The number 
of galaxies with Lya EW > 20A is 87 and this number drops 
to 27 for galaxies with EW > 50A. Such a distribution is par¬ 
tially due to the fact that our sample does not select most of 
the galaxies with very high EW (EW > 50A), whose contin¬ 
uum is usually very faint (~27 mag) and which are usually de¬ 
tected in NB imaging/emission-line selected surveys. This se¬ 
lection effect means that we are probing Lya emitting galaxies 
with brighter continuum compared to galaxies selected based 
on NB/emission-line surveys. We also show the distribution of 
Lya luminosities, L(Lya), for SFGl and LAEs in the middle 
panel of Figure |6] The median L(Lya) of the SFGl sample 
is ~10^® ergs/s, while the median L(Lya) of the LAE sam¬ 
ple is ~10^’ ergs/s. The LAEs discovered in the VUDS sur¬ 
vey are significantly less luminous than the Lya emitters found 
by most N B/emission-line surveys (-^10^^^ ergs/ s) at similar red- 
shifts (e.g.. lGuaita etal.l2010HHagen et al .120141) . The fraction of 
LAEs as a function of redshift is shown in the right panel of Fig¬ 
ure |6] Approximately 10% of SFGs at these redshifts are strong 
Lyg emitters, con sistent with previous Lya studies at z ^ 2 (e.g., 
iReddv et al.ll2008n . and in agreement with the general scenario 
that the number of LA Es increases as redshift increases reaching 
-30-40% atz^6 (e .g., [^k etal.l2nirtrCurtis-Eake et^UOlit 
ICassata et al.ll201^ . 

The composite spectra of the SFGn, SFGl, and LAE sam¬ 
ples are shown in Figure |7] They are generated by stacking the 
rest-frame shifted individual spectra which are rebinned to dis¬ 
persion of 2A per pixel (2A - 5 A/( 1 H-z), where -5 A is the native 
pixel scale of the LR Blue grism in the observed frame), and then 
the entire spectrum is normalized to the mean flux in the wave¬ 
length range of 1350-1450A before taking the median value of 
all galaxies in the composite at each wavelength. Figure|2]shows 
the most prominent emission and absorption lines, whose de¬ 
tailed analysis will be conducted in a future paper. Here, it is 
worth pointing out that the median UV spectral slope (between 
1400-2400A) does not show any appreciable difference among 
these three sub-samples as shown in the bottom panel of Fig¬ 
ure |7] 

The galaxies with strong Lya in emission, usually selected 
using NB imaging or serendipitous spectroscopy, show different 
physical properties compared to galaxies with no Lyg in emis¬ 
sion or LBGs (e.g iGawiser etalj|200d 120071: iLai et al.l 120081: 
ICassata et al.ll201 ih . These two galaxy populations (LAEs and 


LBGs) are selected differently so the differences seen when 
comparing the two populations are mainly due to selection ef¬ 
fects. Here, we use a UV-selected sample (selected based on 
UV continuum/magnitude) to investigate stellar population dif¬ 
ferences between Lya emitting galaxies and non-Lya emitting 
galaxies. A common selection for Lya emitters as well as non¬ 
emitters could reveal real intrinsic differences between these 
two classes of objects. Figure |8] shows the comparison between 
SFGn, SFGl and LAEs. The total number of galaxies identified 
as SFGl and SFGn is large enough (457 SFGn, 397 SFGl, 87 
LAEs) to get robust statistics. On average, SFGl (and LAEs) 
are less dusty, and have lower SLR compared to SFGn. These 
differences are small, but they are significant because of the 
large sample used in this study. To understand the significance 
of these differences we use KS-statistics on these distributions. 
The Pks ~ 10 value for Es(B-V) and SLR distributions im¬ 
plies that the null hypothesis (similar distributions) is rejected 
at >99.9% or >3cr level. The differences in stellar mass, SSFR, 
Mi 5 oo, and /3spec for SFGn and SFGl/LAEs are not statistically 
significant (i.e., the null hypothesis is rejected only at <2.5cr 
level). It is important to note that the KS-test does not consider 
errors/uncertainties in two quantities while comparing their dis¬ 
tributions. We tested the robustness of these results by construct¬ 
ing 100 bootstrap samples by randomly extracting values from 
the confidence intervals of the best-fit stellar parameters. We re¬ 
ran the KS test on each of these 100 artificial samples with a 
goal to compare the stellar population parameters of the LAE 
and non-LAE populations. For Es(B-V) and SLR, a null hypoth¬ 
esis was consistently ruled out at ~3cr level, while for other pa¬ 
rameters, the null hypothesis rejection is at much lower signifi¬ 
cance (<2cr). We also compare the median values of the best-fit 
stellar parameters for LAEs and non-LAEs, and we find similar 
significance for the median values of Es(B-V), SLR, and other 
parameters in these bootstrap samples. These artificial samples 
show that even when accounting for the uncertainty in best-fit 
parameters, a significant difference exists for Es(B-V) and SLR, 
such that galaxies with Lya emission tend to be less dusty, and 
lower in SLR than galaxies with weak or no Lya emission. 

Table [1] summarizes average spectral and photometric prop¬ 
erties of our SFGn, SFGl, and LAE samples. Our results show 
that SFGl (and LAEs) have — on average — lower dust content 
than SFGn. The differences in these two samples is statistically 
significant based on the K-S test results (Pks ~ 10“^), which says 
that distributions of Es(B-V) are different for SFGn and SFGl. 
Though the difference in the distributions of these populations is 
statistically significant, the difference between their median val¬ 
ues is small as shown in Table[T] The median Es(B-V) values for 
SFGl or LAEs and SFGn differ only by -0.05 or -0.1 (i.e., by 
a factor of -1.3 or 2) which is smaller than or similar to the la- 
scatter in these distributions (-0.1). This small difference in the 
dust content of SFGn and SFGl or LAE galaxies is consistent, 
within the large scatter, with the small difference in the median 
values of ySspec, and a small difference in //phot values (-0.1-0.2) 
between these populations. Such a small difference between the 
Es(B-V) values for SFGl or LAEs and SFGn is also consisten t 
with the observations at higher redshifts. [Pentericci et al.l (l2010|) 
found small differences between Es(B-V) values for their sam- 
ples of LBGw and LB Gr at Z'-3. A similar study at z-4 by 
iPentericci et al.l (l2007h agrees with this Es(B-V) trend for galax¬ 
ies with and without Lya in emission. 

Figure |8] also shows that SFGl (and LAEs) have — on aver¬ 
age — lower SED-based SFRs than SFGn. The KS test results 
show that distributions of the SFGl and SFGn samples are sig¬ 
nificantly different i.e., Pks ~ 10 which is as significant as the 


Article number, page 11 ofll9l 


























A&A proofs: manuscript no. ms 



Fig. 6. (Left) Rest-frame EW distribution for SFGs with Lya in emission (positive EWs) and Lya in absorption (negative EWs). The vertical 
dashed line indicates EW = 20A, the adopted lower limit in Lyor EW for LAEs (or strong Lya emitters). The solid curve is the KDE of the 
distribution. (Middle) Distribution of Lya luminosities for the SFGl (black) and the LAE (red) samples. The KDE of the distribution is shown by 
a solid curve. (Right) Fraction of LAEs as a functi on of redshift for our sample. Approx imately 10% of SFGs at these redshifts are strong Lya 
emitters, consistent with previous Lya studies (e.g.. lReddv etal.l[2008l : lGassata etalJPOl.'ih at similar redshifts. 



1200 1400 1600 1800 ,2000 2200 2400 
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Fig. 7. Stacked (median-combined) rest-frame UV spectra for SFGn, SFGl, and LAEs. The strong absorption and emission features are marked by 
the vertical dotted lines. The SEGn spectrum shows strong absorption lines (Si, C, O, Ee, Al), but these absorption lines get weaker with increasing 
strength of Lya-emission. The bottom panel shows similarity in the spectroscopic UV slope for these samples. 


Es(B-V) comparison. The median SFR for the SFGl or LAF 
sample from the SFD-fitting technique is 10* ~ 18 Mq yr ' or 

10 * *^ ~ 15 Mq yr *, while the median SFR for the SFGn sam¬ 
ple is lO*'^’ ~ 23 Mq yr“*. Therefore, the median SFR values for 
SFGl or LAFs and SFGn differ only by a factor of ~1.3 or 1.5, 
which is much smaller in comparison to the Icr scatter of these 
distributions (~0.4 dex or a factor of ~2). To assess the effect o f 
SFD-fitting parameters on the SFRs (e.g.. lKusakabe et^l2015h . 
we also compute SFR using the UV lumino sity (Misnn) correcte d 
for the dust using theyS-Aigoo relation from iMeurer et al.l(Il999l) . 
The median values for the UV-based SFRs (SFRuv) are shown 
in Table [T] and they are very similar to the SFD-based SFRs. 


The median SFRuv for the SFGl or LAF sample is 10* *® ~ 15 
Mq yr“* or 10 M® yr“*, and for the SFGn sample is 

10 * 23~ 17 M® yr"*. Therefore, the difference between SFRuv 
values of these two populations is very small and agrees very 
well with the difference quoted for the SFD-based SFRs. This is 
consistent with the fact that the UV absolute magnitudes (M 1500 ) 
of galaxies with and without Lya in emission are very similar 
as shown in the Table [T| This result is consistent with similar 
studies at z ~ 3-4 (e.g.. lPentericci et^l2007il2010h . 

Figure [8] does not show any significant difference between 
the median values of stellar mass (~0.1 dex or a factor of ~1.2). 
This result is in contrast to higher redshift studies at z ~ 3^ (e g-, 
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Fig. 8. Comparison between SFGn (black histograms), SFGl (red histograms), and LAEs (blue histograms) as a function of stellar parameters. 
The total number of galaxies identified as SFGn, SFGl, and LAEs are shown in the legend. The median value of the distribution is indicated by a 
vertical dashed line. The solid curves are KDE of the distribution. The grey shaded regions cover ± Icr dispersion/scatter, which corresponds to the 
range between the 16th and the 84th percentile values, of the SFGn distribution. The statistical significance based on the KS test (red for SFGl, 
and blue for LAEs with respect to SEGn) is shown to confirm any correlation. 


iPentericci et al.ll2007l 1201 Oh . where they find that galaxies with¬ 
out Lya in emission are more massive (by a factor of ~2-5) com¬ 
pared to galaxies with Lya in emission. The difference between 
these studies could be due to the sample selection, in that the 
stellar mass probed in those studies by the galaxies at z ~ 3^ are 
much lower (down to ~10* Mq) compared to the VUDS sam¬ 
ple (down to ~10® Mq). The lowest mass regime at z~3-4 is 
mostly populated by galaxies with Lya in emission which results 
in this large difference in stellar masses for Lya emitters and 
non-emitters at these redshifts. We do not probe this low mass 
regime between 10* and 10® Mq at 2<z<2.5 because of our 
continuum magnitude limit of Ub ^ 25 mag. Therefore, we con¬ 
clude that, within the luminosities probed by VUDS, we do not 
see any significant difference between the stellar mass of galax¬ 
ies with and without Lya in emission. We will investigate any 
redshift evolution in these correlation in our future VUDS stud¬ 
ies, which will extend these measurements to higher redshifts 
(z>3). 

5. Lyor EW and Stellar Population Properties 

To further investigate properties of SFGn and SFGl, we ex¬ 
plore the correlation between the rest-frame Lya EW and stel¬ 
lar population parameters, such as Misoo, UV slope, Es(B-V), 
SER, stellar mass, and SSER. Eigure|9] shows the correlation be¬ 
tween Lya EW and best-fit SED/spectral parameters. It should 
be noted that the Spearman correlation is measured from all 


galaxies as well as from the median-binned values. The Es(B-V) 
and SER show moderate-to-weak monotonic correlation (p — 
0.2 and Psc < 10 "^) with Lya EW. The significance of p as mea¬ 
sured by the number of standard deviations by which the sum- 
squared difference of ranks deviates from its null-hypothesis ex¬ 
pected value is ~5cr for all Es(B-V) values and ~6cr for all SER 
values. The stellar mass, Mijoo, and /? show much weaker corre¬ 
lation with EW (p < 0.1). The significance of p as defined above 
is <3cr for all values of these parameters. The trends observed 
between the best-fit stellar parameters and EWs are in general 
agreement with the median values of the stellar parameters ob¬ 
tained from the distributions of SEGn and SEGl in Eigure[8] 

With increasing rest-frame Lya EW, we see a weak but 
significant trend in that, galaxies are less dusty, and less star¬ 
forming, as indicated by a decrease in SER by about ~0.2 dex 
from non-LAEs to LAEs. While the average EW varies from 
-38.OA to 120A, the median Es(B-V) decreases from 0.20 to 
0.15, the median SER^^o decreases from 28 Mq yr“' to 16 Mq 
yr“', the median stellar mass decreases from 6.2x10 ® Mq to 
3.9x10® Mq, and the median SSER decreases from 5.9x10 ® 
yr ' to 6.3x10 ® yr '. All these differences in stellar parame¬ 
ters as a function of Lya EW are small compared to the scatter 
observed in each EW bin, as measured by Icr dispersion (see 
Eigure|9]l. 

These results demonstrate that there is a statistically signif¬ 
icant correlation between Lya EW and stellar parameters, such 
as SER and Es(B-V), while no significant correlations are found 
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Table 1. Average properties of SFGn, SFGl, and LAEs derived from the SED-fitting 


Parameters'* 

SEGn'’ 

SEGl** 

LAEs*' 

Ngal 

457 

397 

87 

Ml500 (mag) 

^0 (0.50) 

-ZU.Z/_Q Q 3 

70 1 (0.43) 

z,u.iO_oo4 (0.57) 

7 ^ +0.07 (0.54) 

ZU. 10 _o 10(0.78) 

M 2300 (mag) 

OA c^+0.03 (0.45) 
z,u.>7U_o 03 (0.55) 

7jo ^c+0.03 (0.36) 
z,u.hj_oo4 (0.55) 

70 (0.38) 

^^■^'^-0.09 (0.66) 

Lya Luminosity (Log [erg/s]) 

— 

40 ^7+0.03(0.52) 
HU.JJ_o.05 (0.75) 

41 AQ+0.06 (0.41) 
^^•^^-0.04 (0.31) 

UV slope (spectroscopic) 

. o-r+0.03(0.56) 
^■■^^-0.03 (0.53) 

0^+0.04(0.66) 
^■■^^-0.03 (0.52) 

, ..+0.10(0.80) 
^•^'^-0.07 (0.58) 

UV slope (photometric) 

. 0^+0.03(0.45) 
^■■^^-0.02 (0.35) 

. 0^+0.04(0.59) 
^■■^^-0.03 (0.41) 

, .^+0.11(0.85) 
^■■^°-0.04 (0.30) 

E,(B-V) (mag) 

^ ^(^+0.00 (0.05) 
^■^^-0.01 (0.10) 

^ ,.+0.01 (0.10) 
^•^'^-0.01 (0.10) 

(X ,(x+0.02(0.15) 
-0.01 (0.05) 

SERsed (Log [Mo -yr^i]) 

. 0^+0.02(0.39) 
-0.02 (0.34) 

, ^^+0.02(0.36) 
^■^0-0 02 (0.34) 

, , q+0.05 (0.40) 
^■^^-0.04 (0.28) 

SERuv (Log [Mo -yr”']) 

. ^0+0.03(0.56) 
^■^■^-0.03 (0.57) 

, ,q+0.04(0.59) 
^•^^-0.04 (0.60) 

0 (0-91) 

^•^^-0.09 (0.65) 

Mass (Log [Mo]) 

Q zs+0-03 (0.47) 
^•^'^-0.02 (0.34) 

Q . 7 + 0.03 (0.46) 
-0.02 (0.36) 

Q c^+0.08 (0.62) 

-0.04 (0.33) 

SSER (Log [yr-'D 

0 07+0.02(0.41) 
-o.:)Z_QQ2 (0.36) 

0 ^ 7 + 0.03 (0.46) 
-O.^Z_q o3 (0.44) 

rv .,+0.06(0.43) 
^•^^-0.08 (0.59) 


Notes. Median values of SED parameters and their quoted uncertainties (1.253 x cr/ -y/N^) are shown. The ±lcr dispersion, corresponding to 
the range between the 16th and the 84th percentile values, are also shown in the parenthesis 
® Galaxies with no Lya in emission (EW < OA) 

Galaxies with Lya in emission (EW > OA) 

Galaxies with Lya EW > 20A 


for stellar mass, M 1500 , and p. The differences we observe are 
small compared to the large scatter in their distributions. This 
outcome is true whether or not we focus on all galaxies with Lya 
in emission or only strong Lya emitters (LAEs). This is in con¬ 
trast wit h results from NB or emis s ion-line selected LAE stud¬ 
ies (e . g ., Gawiseret^ 20061 l2007t iLai et al.ll2008l : iGuaita et^ 
I 2 OI It Etagen etal.ll2014l) . These studies find lower dust, lower 
stellar mass and higher SSER for galaxies with strong Lya in 
emission. The LAE samples used in these emission-line/NB- 
selected surveys have Lya luminosities in the range of ~10"^^ 
ergs/s, which is an order of magnitude more luminous than 
VUDS LAEs, as shown in Eigure | 6 ] It is also worth emphasiz¬ 
ing that VUDS is a UV continuum selected survey and galaxies 
are targeted based on their photometric redshifts and not on the 
strength of the Lya line. We have 27 galaxies beyond rest-frame 
Lya EW = 50A, so by selection, our sample does not include 
extremely low mass (~10^ Mq) galaxies usually selected from 
emission-line/NB techniques. Moreover, among the emission- 
line/NB-selected LAEs, those that have lower Lya luminosi¬ 
ties, have similar SER as normal star-forming galaxies but lie 
on the lower mass en d of the star-forming main sequence (e.g., 
iKusakabe et al.ll^l5h . Hence, the difference in the sample se¬ 
lection could play a significant role in such a comparison as 
we could be comparing two intrinsically different populations 
of galaxies. 


The comparison of stellar population properties of galax¬ 
ies with and without Lya emission has also been done on 
UV-selected galaxies at z ^ 2 and high er, and our res u lts are 
broadly consistent with the literature. iPentericci et al.l (1201 Ol) 
found weak/no correlation between SERs of galaxies with and 
without Lya emission at z ^ 3 but we find a correlation (with a 


small difference in SER) possibly because of the large statistics 
in the VUDS sample. With a smaller number of galaxies (130) in 
the Pentericci et al. sample such a detection was difficult to ob¬ 
serve. Our finding of a very small decrease in the SED-based 
dust content (i.e., ETB-V)) of th ese two popul ations agrees 
well w ith the IPentericci et alJ (I 2 OIOI) study. At z - 3. lKornei et al.l 
(I 2 OIOI) find that LAEs have lower Es(B-V) and SERs compared 
to non-LAEs. The difference in these two stellar population pa¬ 
rameters for LAEs and non-LAEs, as found by Kornei et al., is 
much larger than what we observe, but show similar decreas¬ 
ing trends in these two quantities. This trend of lower dust and 
SER for strong Ly a emitters ha s also been observe d at z - 2 by 
iReddv et"^ (1200^ . The study of iKornei et alJ (1201 Ol) agrees well 
with our result that there is no significant difference between 
stellar masses of LAEs and non-LAEs. Therefore, based on the 
comparison between SEGl/LAEs and SEGn, as well as Lya EW 
correlations with stellar parameters, we suggest that our LAE 
sample shows a small decrease in Es(B-V) and SER but oth¬ 
erwise their stellar populations are not very different from the 
galaxies that do not show Lya in emission. 

iLai et al.l (l2008l) studied NB-selected LAEs at z - 3 by divid¬ 
ing them into Spitzer/IRAC-detected (brighter than m 3 6 = 25.2 
mag) and IRAC-undetected (fainter than m 3 6 = 25.2 AB mag) 
samples. They found that IRAC-undetected LAEs were less mas¬ 
sive and younger compared to IRAC-detected LAEs, concluding 
that IRAC-detected LAEs are a more evolve d population similar 
to bri ghter/massive UV-selected LAEs (e.g.. lShaplev etal]l200lL 
l2003h . To investigate this trend for VUDS LAEs, we examined 
our LAEs as a function of magnitude in the Spitzer/IRAC 3.6jum 
band. We find that 74% of SEGn, 66 % of SEGl, and 53% of 
LAEs have detections in the Spitzer/IRAC 3.6jum band down to 
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Fig. 9. Correlation between rest-frame Lya EW and best-fit stellar population parameters. The vertical dashed line indicates EW = 20A, the adopted 
lower limit in Lyn equivalent width for LAEs. The size of the last bin is increased with respect to that of the other bins to increase the number of 
galaxies in that bin. The light grey points are individual measurements while the hlack circles are binned medians. The error bars in x illustrate 
the sizes of the bins, while the errors in y are ±1 <t scatter (dashed error bars) corresponding to the range between 16th and the 84th percentile 
values within each bin, while smaller solid error bars are the errors on the median values (1.253 x <t/ ^yNgai). The statistical significance from 
the Spearman correlation coefficient for all galaxies (median-binned values) is shown to confirm these correlations. Average uncertainties in these 
measurements are plotted in the bottom-right comer. The median uncertainties in the EW measurement varies as a function of EW as shown by 
different x error bars. 


~25 mag. The LAE fraction with IRAC dete ction in our sam - 
ple is much larger than the ~30% found in iLai et al.l (l2008h . 
which is a direct consequence of the difference in the sample 
selection between the two studies. The NB-selected LAEs typ¬ 
ically have very faint continuum in all optical and NIR bands 
as they are mostly lower mass/younger galaxies compared to 
normal star-forming galaxies. IRAC-detection is only possible 
in NB-LAEs if they host more evolved population which makes 
continuum brighter in the red bands. We do not expect many NB- 
selected LAEs to be more evolved/massive and brighter in IRAC 
bands, therefore their IRAC-detected fractions are more likely 
to be small c ompared to other s tar-forming galaxies. A similar 
study at z ^ 2 (iGuaita etal.ll201 ih . also found that ~25% of NB- 
selected LAEs had IRAC detection down to m 3 g ~ 24.5 AB mag. 
We find a higher detection rate (51 %) of IRAC-detected LAEs 
in our sample using the same limiting magnitude of 24.5 mag as 
Guaita et al. EigurefTOlshows the comparison of IRAC-detected 
and IRAC-undetected galaxies in all three samples as a function 
of stellar mass, and SER. The stellar mass correlations show that, 
in all three samples, IRAC-detected galaxies are more massive 
than IRAC-undetected galaxies. These correlations are statisti¬ 
cally significant based on their Pks values < 10“^. We do not see 
any significant correlation for SER in all three samples, both in 
terms of Pks and their median values. 

These comparisons show that the VUDS LAE sample has 
a larger fraction of galaxies with IRAC detection compared to 
NB-selected LAEs at z ^ 2-3, and that IRAC-detected LAEs are 


more likely to host evolved populations compared to IRAC- 
undetected LAEs. A large fraction of UV-selected LAEs have 
similar masses as SFGn galaxies but with lower SERs suggesting 
that these LAEs host more evolved populations. The low frac¬ 
tion of IRAC detection in NB-selected LAEs implies that these 
galaxies are less evolved compared to UV-selected LAEs. 


6. Results and Discussion 

We use a large spectroscopic sample of 854 SEGs at 2 <z<2.5 
from VUDS to investigate their spectral and photometric prop¬ 
erties. The VUDS spectra were used to measure the UV spectral 
slope and Lya EW. The SED fitting on multi-wavelength pho¬ 
tometric data using Le PHARE provided best-fitting stellar pa¬ 
rameters, such as stellar mass, M1500, Es(B-V), SER, and SSER. 
The VUDS observations extend to a fairly large area of ~ 1 deg^, 
which implies that our sample spans a large range in SER (~3- 
150 Mo yr-') and stellar mass (-5x10*^ to 10“ Mo). The VUDS 
2<z<2.5 sample covers a substantial range in UV absolute 
magnitude around M* (+1 mag), the characteristic magnitude 
at these redshifts. We divide SEGs into three sub-groups, SEGn, 
SEGl and LAEs, based on their rest-frame Lya EW. We find 
that the fraction of LAEs is ~10% a t these redshifts, w hich is 
consistent with previous studies (e.g.. lReddv et al.l[2008h . and is 
in accordance with the gen eral picture where the Lya fraction 
increases with redshift (e .g., IStark etani2010l : ICurtis-Lake et al.l 
l 2012 l : ICassata etani2015h . The Lya fraction apparently starts to 
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Fig. 10. Comparison between Spitzer/IRAC 3.6/rm-detecteci (red; brighter than ~ 25 mag) and IRAC-undetected (blue; fainter than m 3.6 ~ 25 
mag) galaxies as function of stellar mass and SFR for SFGn, SFGl, and LAEs. The fraction of galaxies detected/un-detected in IRAC are shown 
in the legend. The solid curve is the KDE of the distribution. The statistical significance based on the KS test is shown to confirm the observed 
correlation. 


decrease at higher redshift (z>6.5; e.g.. lPentericci etani2011h . 
where large amounts of neutral hydrogen start to affect the visi¬ 
bility of Lya line indicating the onset of the reionization epoch. 

The UV spectral slope /? can be used t o derive the dust at¬ 
tenuation for local s tarburst galaxies (e.g.. ICalzetti etal1ll994[ 
iLeitherer et al.l [199^ and has bee n adopted as a dust indica¬ 
tor for high redshifts galaxies (e.g.. [Noll et al] |2004 ; iHathi et all 
l2008bl; [pinkelstein et al.l l20T^ lBou wens et alJl2014 ). The main 
reason is that the intrinsic UV slope depends only weakly on 
metallicity and stellar populations for star-f orming galaxies (e.g., 
iHeckman et al.|[T9^ILeitherer et al J 1999t) . In addition,/? shows 
a strong correl ation with the LTR-to-Lriv ratio, a sta ndard dust 
indicator (e.g.. iMeurer et al.lll999l: iReddv et al.ll201(lh . On aver¬ 
age, the spectroscopic UV slope for SFGs at 2<z<2.5, mea¬ 
sured using VUDS spectra, is comparable to the photometric /? 
measured using multiple photometric bands, and have smaller 
measurement uncertainties. The measured UV slope — spectro¬ 
scopic and photometric — is in accordance with the evolutionary 
trend of fi with redshift, which shows that lower redshift galaxies 
have, on aver age, redder UV slope compared to high er redshift 
galaxies (e.g., iBouwens etalJl201^[Hathi etalj|20ll) . With dust 
being the major factor influencing the UV slope, this implies that 
lower redshift galaxies have more dust compared to higher red¬ 
shift galaxies. We use the spectroscopic p measurements to ex¬ 
plore its correlation with Misoo and stellar mass. Comparing/! to 
the stellar mass we find a significant correlation, in the sense that 
massive galaxies are redder, which is in general agreement with 
the higher redshift measurements of iFinkelstein etal](l2012h . We 


find no correlation between M 1500 and /?. This result is at vari- 
ance with what found by some authors at higher redshifts (e.g., 
iBouwens etal.l20i^ [Rogers et al.l20'r4l) . but it is similar to other 
studies which do not find any significant correlation between P 
and Mi-snn Finkelstein et al.l 120121 : ICaste llano et al.l 1201^ 

iFteinis etal.ll2013h nr is vital to note that this correlation could 
be affected by the biases and/or differences in the way differ¬ 
ent authors measure M 1500 and/or p but most importantly the 
dynamic range in M 1500 plays a crucial role. Extensive studies 
exploring a larger range in M 1500 and stellar masses and investi¬ 
gating different biases in these correlations will shed more light 
on the true nature of these relations and provide a better physical 
understanding. 

There are several studies investig ating the stellar population 


ofUV-selected LAEs at 2-3 (e . g..lStolev etalJl2bo 11120031 
lErb etalJl2()()5 iReddv et all 120081: iKOTnei etal.ll201oir We find 


that LAEs in our sample have lower dust content, lower SER 
and lower SSER compared to non-emitters. These differences 
are very small compared to the large scatter in these SED-based 
paramet ers but they are sig nific ant because of the large sam¬ 
ple size. IShaplev et al.| (1200 ll) and iKornei et al.l (1201 Oh found that 
LAEs at z - 3 have lower SERs, lower dust and lower SSER com¬ 
pared to LBGs, which is consistent with our observed trends, 
but we find much smaller differences in these stellar parameters 
for our galaxies. We find only weak or no correlation of stellar 
mass with Lya EW for a sa mple with median L og(Mt) = 9.64 
and median Log(SER) = L32. lKornei etalJ(l2010h studied a large 
(~300) sample of UV-selected galaxies, with and without Lya 
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emission, with a median stellar mass of Log(Mt) = 9.92 and a 
median SFR of Log(SFR) = 1.57, both slightly higher than for 
our sample. They found that there is no significant difference in 
the stellar mass between these two populations in accordance 
with our weak-to-no correlation of stellar mass with the Lycr 
EW. The results obtained here show no strong correlation or 
large differences between stellar parameters as a function of the 
Lyg EW, which is consistent with the small sample study of 
[Reddy etaP (l2008h using UV-selected galaxies at z-2-3. It is 
clear that the properties of galaxies with and without Lya in 
emission are pretty similar suggesting that the two populations, 
at least for typical (~L*) galaxies and to the level of detail we 
are able t o prob e, are roughly comprised of similar galaxies. 
lErb et al.l (l2006h used composite spectra of UV-selected LAEs 
and non-LAEs at z ^ 2 and found that galaxies with strong Lya 
EW had lower stellar masses. Our results on SER and dust are 
broadly consistent with previous studies, but the lack of a strong 
correlation between stellar mass and EW needs to be explored 
with a sample that spans a larger dynamic range in stellar mass 
as lower mass galaxies could strongly impact the stellar mass 
correlation between Lyo- emitters and non-emitters. 

In comparison to UV-selected LAEs, LAEs selected based 
on the emission-line/NB technique at z^2 and beyond have 
very low masses ~10^ Mp compared to Mp for non- 

LAEs (e.g., Gawi^er_e^I.|20^[Fiifel^ste^ et aI.l2007HLai et aP 

1 IVargas el 


l2008t iGuaita et all 201 it Ivargas et al.ll20I4 ). Also, these LAEs 
are much more luminous in Lya (~10'*^ ergs/s) than VUDS 
LAEs, which are on average an order of magnitude less lu¬ 
minous (~10^* ergs/s). We do not find such low mass, high 
Lya luminosity LAEs in our sample mainly because these 
strong emitters have low continuum luminosities and are not se¬ 
lected by the VUDS magnitude selection. In fact, such galaxies 
woul d be missed in almost all UV continuum selected sample s 

te.g.- l^anlev et af 1I200II: iReddv et al.ll2008l: iKornei etal.ll201(]h . 

as emission-line/NB LAEs probe a different lumi nosity range 
compa red to magnitude-limited samples. However. iKornei et al.l 
(I 2 OIOI) have argued that when NB-selected LAEs are restricted to 
luminosities probed by continuum selected samples, both sam¬ 
ples are statistically very similar, which is consistent with the re¬ 
sult put forward by the study of IVerhamme et all (l2008h . where 
bright NB-selected LAEs represent the same population as con¬ 
tinuum selected LAEs. Therefore, UV-selected LAEs and NB- 
selected LAEs need to probe similar luminosities for proper 
comparison between these two populations. 

In the literature there are two possible scenarios under which 
galaxies emit Lya during the evolution process . Studi e s of faint 
NB-s elected LAEs (e.g., iGawiser et al'l 120061 l2007l iLai et al.l 
I 2 OO 8 I) have argued that LAEs represent the beginning of an evo¬ 
lutionary sequence of galaxy formation. According to these stud¬ 
ies, the main reason behind this argument is that LAEs have 
lower dust content, lower stellar mass, younger stellar ages, 
and higher SSERs compared to non-LAEs which means that 
LAEs are building up their stellar mass at a faster rate than 
non-LAEs through mergers and/or star-formation episodes. On 
the o ther hand, studies of brighter continuum sele cted LAEs 
(e.g., IShanlev et al.l 120011 l2003t IKornei et^ l2010h have sug¬ 
gested that LAEs represent a later stage in the evolutionary se¬ 
quence. These studies argue that LAEs have lower dust content, 
lower SFR, older stellar ages, and lower SSFRs compared to 
non-LAEs which means that LAEs are more quiescent than non- 
LAEs. This could be the result of strong outflows from super¬ 
novae and massive star winds, which expel both gas and dust 
from a young, dusty non-LAE. These arguments raise a ques¬ 
tion: are these two scenarios two different phases of the evolu¬ 


tion proc e ss or d o these two possibilities describe a single phase? 
iLai et al.l (l2008l) investigated the stellar populations of 162 NB- 
selected LAEs at z ^ 3.1 by dividing them in two groups based 
on their Spitzer/IRAC flux. They find that ~70% of the LAEs are 
undetected in 3.6/zm down to 25.2 mag. Based on their stacking 
analysis, they find a clear difference between these two divided 
samples. The average stellar population of the IRAC-undetected 
sample had an age of ~200 Myr and a mass of ~3xl0^ Mq, 
consistent with the scenario that LAEs are mostly young and 
low-mass galaxies. On the other hand, the IRAC-detected LAEs 
were on average significantly older and more massive, with an 
average age of ~1 Gyr and mass of ~10'° Mf7). The s tellar pop- 
ulations of the IRAC-detected sample of iLai et^ (l2008li are 
very similar to those of continuum selected bri ghter LAEs (e.g., 
IShartlev et al.l 1200 iL 120031 : IKornei etal.l l2010t) . Similar results 
are also found in the study of Guaitaetal.1 (I 2 OIII) at z-2.1, 
in that, the fraction of IRAC-detected LAEs is ~25% down 
to 24.5 mag and these galaxies are more massive/evolved than 
those that are faint in IRAC. We find larger fraction of IRAC- 
detected LAEs (~50%) compared to NB-selected LAEs (~30%) 
at z - 2-3, which could imply that UV-selected LAEs are more 
evolved compared to NB/emission-line selected LAEs. Also, the 
IRAC-detected LAEs have higher stellar mass, similar to IRAC- 
detected non-LAEs, but with lower SFR which is consistent with 
the more evolved nature of these galaxies. The scenario in which 
LAEs are more evolved than non-LAEs is plausible if strong 
outflows in non-LAEs destroy or remove dust and gas from the 
galaxies, allowing Lya photons to escape. Our results from di¬ 
rect comparison of galaxies with and without Lya emission show 
small differences in stellar populations, such that galaxies with 
Lya emission have lower SFR and dust content than non-Lya 
emitting galaxies. This result broadly agrees with previous stud¬ 
ies involving UV-selecte d galaxies at z - 3 (e.g.. IShanlev et all 
l2003t IKornei et al.l l2010li . though we observe small difference 
between physical paramete rs of Lya emitters and non-emitters. 
Past studies of LAEs (e.g..lyerhamme et al.ll2008t iNilsson et al.l 
l201lt lAcauaviva et al.l 120121) . suggest that LAEs have a large 
range in stellar population properties implying that a strong L AE 
phase either is a long duration phase (~1 Gvr: lLai et al.ll200^ or 
is recurring in SEGs. 

Studies focusing on comparing stellar populations of galax¬ 
ies with and without Lya emission have shown that various ef¬ 
fects, including the sample selection, continuum and Lya lumi¬ 
nosities, measurements of stellar parameters, and SED fitting 
assumptions play an important role in how we compare Lya 
emitters and non-emitters. It is also essential to understand that 
these comparisons and correlations are not universal and could 
change with redshift. Extending these studies to higher redshifts, 
where the fraction of Lya emitters increases, is essential to gen¬ 
erate a self-consistent picture of how LAEs and non-LAEs form 
and evolve. In future studies, we will extend such an analysis 
to z - 3-6 using the VUDS data, to better understand how these 
LAE properties/trends evolve with redshift from z ~ 6 to z ~ 2. 

7. Summary 

In this paper, we have investigated the spectro-photometric prop¬ 
erties of a large sample of SFGs at 2<z<2.5 that were se¬ 
lected from the VIMOS Ultra-Deep Survey in the ECDES, COS¬ 
MOS, and VVDS fields. These galaxies were targeted because 
of their photometric redshifts, and are therefore UV continuum- 
selected galaxies. The VUDS spectra were used to measure the 
UV spectral slope (J3) and Lya EWs, while we used deep multi¬ 
wavelength observations in these extensively observed fields to 
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derive physical parameters (stellar mass, SFR, Es(B-V), M 1500 , 
SSFR) from the SED fitting process. We compared Lya emitters 
and non-emitters using these parameters and also explored cor¬ 
relations of these parameters with Lya EW. Our results can be 
summarized as follows: 

• We obtain reliable measurements of spectroscopic UV 
slopes, which are — on average — similar to the photometric 
measurements of (3, and have smaller measurement uncertain¬ 
ties. The median values of jSspec and ySphot are consistent with 
each other and with the general picture in which UV slopes get 
redder with decreasing redshift, implying higher dust content at 
lower redshifts. We find a significant correlation between (3 and 
SED-based dust indicator Es(B-V). We observe no correlation 
between (3 and M 1500 , while a strong correlation between (3 and 
stellar mass is observed. These results are consistent with higher 
redshift observations. 

• For a proper comparison, we divide these SFGs into 
three subgroups based on their EWs: SFGn (EW <0A), SFGl 
(EW > OA), and LAEs (EW > 20A). The LAEs make up -10% 
of the total SFG sample at 2<z<2.5, which is consistent with 
previous observations. 

• We find that at 2 < z < 2.5, within the luminosities probed, 
the SFGl (and LAE) sample has slightly lower Es(B-V) and 
SFRs compared to the SFGn sample. These differences are small 
but statistically significant. It is important to note that we are able 
to probe these small differences in physical parameters because 
of our large SFG sample. We do not find any or find weaker sig¬ 
nificant differences in stellar mass, M 1500 , and (3 for these two 
samples. We find similar results when we compare Lya EW and 
physical parameters. These results indicate that the properties 
of galaxies with and without Lya in emission are remarkably 
similar suggesting that the two populations, at least for typical 
(~L*) galaxies and to the level of detail we are able to probe, are 
roughly comprised of similar galaxies. 

• When we divide the LAEs based on their Spitzer/IRAC 
3.6/im magnitudes, the fraction of IRAC-detected LAEs (-50%) 
is much greater than the fraction of IRAC-detected NB-selected 
LAEs (-30%) at z^2-3. This could imply that UV-selected 
LAEs host a more evolved stellar population compared to 
NB/emission-line selected LAEs. Based on differing stellar pop¬ 
ulation results for LAEs at various redshifts, we cannot rule out 
multiple or recurring Lya emitting phases for SFGs. 

Future studies of higher redshift (z - 3-6) galaxies from 
VUDS, where the fraction of Lya-emitting galaxies increases 
substantially, will help us to better understand the LAE popu¬ 
lation at these redshifts, as well as to see how these correlation 
evolve with redshift. 
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